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ABSTRACT 

HH 211 is a highly collimated jet originating from a nearby young Class 
protostar. Here is a follow-up study of the jet with our previous observations 
at unprecedented resolution up to ~ 0!'3 in SiO (J = 8 — 7), CO (J = 3 — 2), 
and SO {Nj = 89 — Tg). SiO, CO, and SO can all be a good tracer of the 
HH 211 jet, tracing the internal shocks in the jet. Although the emissions of 
these molecules show roughly the same morphology of the jet, there are detailed 
differences. In particular, the CO emission traces the jet closer to the source than 
the SiO and SO emissions. In addition, in the better resolved internal shocks, 
both the CO and SO emission are seen slightly ahead of the SiO emission. The 
jet is clearly seen on both sides of the source with more than one cycle of wiggle. 
The wiggle is reflection-symmetric about the source and can be reasonably fitted 
by an orbiting source jet model. The best-fit parameters suggest that the source 
itself could be a very low-mass protobinary with a total mass of ~ 60 Mjup and a 
binary separation of ~ 4.6 AU. The abundances of SiO and SO in the gas phase 
are found to be highly enhanced in the jet as compared to the quiescent molecular 
clouds, even close to within 300 AU from the source where the dynamical time 
scale is < 10 yrs. The abundance enhancements of these molecules are closely 
related to the internal shocks. The detected SiO is either the consequence of 
the release of Si-bearing material from dust grains or of its formation via gas 
chemistry in the shocks. The SO, on the other hand, seems to form via gas 
chemistry in the shocks. 
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Introduction 



Jets from protostars represent one of the most intriguing signposts of star formation. 
They are beheved to be launched from accretion disks around the protostars. They are 
highly supersonic, coUimated, and ballistic. Yet, they are often found to show wiggles in 
their trajectories . These wiggles of the jets could be due to a precession of the jets (e.g.. 



Raga et al. 



1993f). or an orbital mot i on of the jet sources (i.e , the protostars that drive the 



jets) (e.g., iFendt fc Zinneckerl Il998t iMasciadri fc Ragal 120021 ). or both. The jet precession 
may be due to a precessi on of the accretion dis ks because of tidal interactions in noncoplanar 
binary systems (see, e.g, iTerquem et al.l 119991 ) and it can give rise to point-symmetric (i.e., 
S-shaped) wiggles, while the orbital motion of the jet sources may arise in binaries an d 
it can give rise to reflect ion-symmetric (i.e., C-shaped) wiggles (IFendt fc Zinneckerl Il998l ). 
Masciadri fc Ragal (|2002[ ) have studied the wiggles of the DG Tau microjet, HH 47, and the 
Serpens radio continuum jet. Similarly, lAnglada et al.l ( 120071 ) have studied the wiggle of the 
HH 30 jet. In both cases, assuming that the wiggles of those jets are due to an orbital motion 
of the jet sources, they have derived orbital parameters and masses that are reasonable for 
pre-main-sequence binaries. For those jets, however, only one side of the jets is seen, and 
thus we can not confirm if the wiggles are really reflection-symmetric. 

This paper is a follow-up study of the HH 211 jet with our previous observations at 
unprecedented resolution u p to ~ ^ ^3 (Le e et al. 2009, hereafter iPaper j ) . obtained with 
the Submillimeter Arrajll] ( Ho et al. 2004). The jet is clearly seen in SiO (J = 8 — 7) on 
both sides of a young Class source with more than one cycle of wiggle, and is thus one 
of the best candidates to study the origin of the wiggle. The wi ggle is re flection-symmetric 
and it seems to be due to an orbital m otion of the jet source (jPaper j ). The wiggle was 
called the C-shaped bending in iPaper a but not here anymore because more than one cycle 
is seen. In order to have a more complete picture of the wiggle, here we also present CO 
( J = 3 — 2) and SO {Nj = 89 — 7%) emissions of the jet in the same observations as already 
described in iPaper 1 By modeling the wiggle, we derive the orbital parameters and then 
discuss the nature of the source. In addition, we also refine the mass-loss rate of the jet with 
the CO emission at high resolution and investigate the possible origins of the SiO, SO, and 
CO emissions in the jet. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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Results 



In the following, we present the morphology and kinematics of the jet in CO (J = 3 — 2) 
and SO {Nj = 89 — 7s) emissions in comparison to those in SiO (J = 8 — 7) emission in 
Paper t . We also derive the density of the jet from the CO emission and then the abundances 
of SO and SiO. The systemic velocity in HH 211 is assumed to b e 9.2 ± .3 km s~^ LSR and 
the distance to the HH 211 is assumed to be 280±40 pc, as in iPaper t . The jet is almost 
in the plane of the sky with an inclination angle of ~ 5°±3°. The eastern component of 
the jet is blueshifted while the western component is redshifted. The eastern component of 
the jet has a P.A. of 116.1°±0.5°, while the western component has a P. A. of 297 .1°±0.5°, 
indicating that the jet has a small-scale bending, ~ 0.5° to the northeast (jPaper J ). 



2.1. Morphology of the jet 

Figure [1] shows the CO and SO maps on top of the Si O map of the jet and the 352 GHz 



continuum map (green contours) of the envelope-disk from lPaper 11 . The CO emission of the 



jet is derived by excluding the CO emission within ~ ± 10 km s ^ from the systemi c, in order 



to avo id the contamination from the shells and internal shells (for more details, see lLee et al. 



2007bl ). The CO emission of the jet is more than a factor of 4 weaker than the SiO emission. 
Although the CO emission shows roughly the same morphology as the SiO emission, there 
are detailed differences (Fig. Wp)- The CO emission shows a chain of paired knots on either 
side of the source SMMl, but with the peaks slightly ahead (i.e., downstream) of the SiO 
peaks for most of the better resolved knots (see, e.g., knots BK2, BK3, and RK2 in Fig. [T]^). 
The CO emission also shows the same curvy structures for knots BKl and RKl, and the 
same reflection-symmetric wiggle of the jet (Figs. [Th, c) as SiO. However, the CO emission 
extends closer in to the source than the SiO emission, with the closest peaks (labeled as 
knots RKO and BKO in Fig. [lb) at ~ O'.'A (~ 100 AU) on either side of the source. The SO 
emission is even weaker than the CO emission, but it also shows knotty structures along the 
jet, like the CO emission (Fig. [1^). Like the CO peaks, the SO peaks also appear slightly 
ahead of the SiO peaks in e.g., knots BK2, BK3, and RK2. However, almost no SO emission 
is seen toward the innermost pair of CO knots. 



2.2. Kinematics along the jet axis 



Figure [2] shows the position- velocity (PV) diagram cut along the jet axis in CO for the 
two inner pairs of knots, e.g., BKO, BKl, RKO and RKl, in comparison to that in SiO. The 
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innermost pair of CO knots, BKO and RKO, allow us to study the jet kinematics closer to the 
source than the innermost pair of SiO knots, BKl and RKl. They are seen associated with a 
broad range of velocities, tracing the internal shocks closer to the source than the S iO knots . 



The four sub-knots of BKl and RKl seen in SiO with a range of velocities (see iPaper 11 ) 
have possible counterparts in CO in the PV diagram (Fig. [2], marked with yellow lines and 
question marks). The separation between the first SiO sub- knots and the innermost pair of 
CO knots is ~ 0''8 (or ~ 220 AU), roughly the same as that in between the sub-knots. This 
suggests that the innermost pair of CO knots and the sub-knots of BKl and RKl are all the 
internal shocks produced b y the sa me mechanism, namely by a periodical variation in the 



jet velocity as suggested in iPaper I . 



The knots BK2 and BK3 are better resolved, allowing us to study the detailed kinematics 
in different emissions in the internal shocks. Figure E] shows the PV diagrams cut along the 
jet axis in CO, SiO, and SO for these knots. In SiO, these knots clearly show a head-tail 
morphology, with the heads associated with a broad range of velocities and thus tra cing 



the internal shocks. With respect to the mean velocity of the jet reported in iPaper 11 . the 



nearside of the heads is blueshifted and the fars ide is redshifted, because o f the sideways 



ejection of the jet material in the shocks (see also ISantiago-Garcia et al.ll2009l ). In the heads 
of these knots, the CO emission appears ahead of the SiO emission and the velocity range of 
the CO emission is smaller than that of the SiO emission. Since SiO emission traces denser 
and warmer regions than CO emission, these differences suggest that the SiO emission traces 
stronger shocks closer to the shock front and the CO emission traces the weaker shocks in 
the downstream. As for the SO emission, only knot BK3 is bright enough for kinematic 
study. In that knot, the SO emission is slightly ahead of the SiO emission with a smaller 
velocity range, tracing weaker shocks than the SiO emission. Also, the SO emission seems 
slightly behind (i.e., upstream) the CO emission and peaks at higher velocity offsets from 
the mean jet velocity than the CO emission, tracing stronger shocks than the CO emission. 



2.3. Column and Volume Densities of the jet 

The column and volume densities of the jet can be estimated from the CO emission of 
the two innermost pairs of knots (i.e., knots BKO, BKl, RKO, and RKl). At the resolution 
of ~ 0''46x0''36, the mean intensity of the CO emission toward these knots is ~ 3.2 Jy 
beam~^ km s~^ (see Fig. [T)d). The abundance of CO relative to H2 is assumed to be 
4 X 10~^ as if the CO gas is formed via gas-phase reactions in an initially atomic jet 



(iGlassgold et al.l Il99ll . see also §4.31) . The kinetic temperature the CO emission can be 



assumed to be 200 K, slightly lower than the kinetic temperature derived from the CO 
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emission in much hig her transitions in the far infrared, which was found to be > 250 K 
( IGiannini et aLll200ll ). Assuming LTE, the excitation temperature of the CO emission can 
be assumed to be the same as the kinetic temperature. Assuming that the CO emission 
is optically thin, the column density of the jet in H2 is ~ 5.5 xlO^*^ cm~^. However, 
the true column density of the jet is actually higher because the jet is spatially unresolved 
perpendicular to the jet axis. With the beam having a size b ~ 0!'45 or 125 AU perpendicular 
to the jet axis and the jet having a deconvolved width < 40 AU (derived from the SiO jet 



m 



Paper t ). the true column density of the jet is A'a ^ A^(125/40) ~ 1.7 x 10^^ cm~^. As a 
result, the volume density of the jet is n > A'a/40AU ~ 2.8 x 10^ cm"^, about two orders of 
magnitude higher than the critical density of CO ( J = 3 — 2) line, which is 3 x 10^ cm~^ at 
200 K. On the other hand, the volume density derived here is one order of magnitude lower 
than the critical density of the SO line, and two orders of magnitude lower than that of the 
SiO line at 200-500 K. Note that the critical densities here are derive d using the coUisio nal 
rate coefficients listed in the Leiden Atomic and Molecular database (jSchoier et al.ll2005l ). 



2.4. SO and SiO abundances in the jet 



The abundances of SO and SiO relative to H2 can be estimated by dividing SO and 
SiO column densities by that of H2 derived from the CO emission, respectively. However, 
since the SO and SiO emissions trace stronger shocks and thus denser material than the CO 
emission, the abundances estimated here are likely to be the upper limits. 

At the resolution of 0''46x0''36, the mean SO intensity toward the innermost pair of SO 
knots (i.e., knots BKl and RKl) is ~ 1 Jy beam~^ km s~^ (see Fig. [T^). The excitation 
temperature of the SO emission is assumed to be the same as that of the CO emission, but 
it could be higher because the SO emission seems to trace stronger shocks than the CO 
emission as seen in knot BK3 in the downstream (see §2.21) . Assuming that the SO emission 
is optically thin, the SO column density is ~ 1.1 x 10^^ cm~^. Therefore, the abundance of 
SO is ~ 2xl0 ~^. It is ~ 400 times t hat in cold quiescent clouds (cores), e.g., TMC-1, which 
is ~ 5 X 10"^ Johishi fc KaifulllQQsh . 



At the resolution of 0''46x0!'36, the SiO intensity toward the innermost pair of SiO 
knots decreases rapidly with the distance from the source from ~ 20 to 5 Jy beam~^ km s~^ 
(Fig. W^ ). Assuming an excitation temperature of 500 K and optically thin emission as in 
Paper 11 . the SiO column density at the SiO knots is (4.4-1.0) x 10^^ cm~^, resulting in a SiO 
abundance of (8-2)xl0~®. Therefore, the SiO abundance is highly enhanced by about five 
orders of magnitude oyer th at in quiescent molecular clouds where the abundance of SiO is 
< 10-11 Jziurvs et al.lll989h . 
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3. An orbiting source jet model for the wiggle 



The reflection-symmet ric wiggle of the jet could be due to an orbital m otion of the jet 
source in a binary system (IFendt &: Zinneckerlll998l : iMasciadri &: Ragall2002l ). If that is the 
case, then by modeling the wiggle of the jet, we could obtain the parameters of the orbital 
motion and then discuss the nature of the source (see §4.ip . To model, we adopt the right- 
handed Cartesian coordinate system as shown in Figure HJ with the +z-axis aligned with 
the eastern component of the jet. For simplicity, the jet source is assumed to have a circular 
orbit in the x-y plane, with the jet propagating ballistically in the 2;-axis at a velocity of 
Vy The orbit has a period Pq, an orbital velocity fo, and a radius Rq. Assuming that the 
orbital motion of the je t source is in the same direction as the rotation motion of the HCO"'" 
envelope-disk (jPaper J ), we have f o > (i.e., positive for right-handed orbital motion). Then 
for a jet lying close to the plane of the sky as in the case of HH 211, the trajectory of the 
jet in the plane of the sky can be approximately given by 



X = -Ji?H(^o-)^cos[0o-^- + tan-i(^^)] 



+r] z 



Ro 



Ro v^' 



(1) 



where 0o is the current phase angle of the jet source in the orbit measured from the +x-axis 
toward the -l-jy -axis dir ection, and rj is to account for the small-scale bending of the jet to 
the northeast (jPaper 11 ) . Deflne the periodic length as 



and the velocity ratio as 



we have 



and thus 



Az = v;Po 



K = 



Ro 



VoPo K,Az 
^tT ~ 271 



(2) 
(3) 
(4) 



X 



+ I cos(0o 1. h tan 



27r 



+77 z 



Az 



Az' 



(5) 



When z ^ tA, have 



271 Z 71 , 



X ~ —KZ COS[(pQ 



Az 2 



(6) 
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which means that the wiggle due to the orbital motion of the jet source has an opening 
angle a ~ 2k. Figure [5] shows the best fit of this model to the SiO, CO, and SO data 
simultaneously. As can be seen, the wiggle of the jet can be reasonably fitted by this model. 
The best-fit parameters are Az = 1530 ± 25 AU (~ 5('46±0:'09), k = 0.0094 ± 0.0007 (or 
0.54°±0.04°), 00 = 155°±13°, and = 0.0098±0.0009 (or 0.56°±0.05°). Thus, Ro ~ 2.3±0.2 
AU, Vn ~ 1.6 ± 0.56 km s~^, and Pq ~ 43 ± 23 yrs, assuming that = 170 ± 60 km s~^ 



as m 



Paper 11 . In this model, the wiggle of the jet has an opening angle of a ~ 1-1°; similar 



to that seen in the observations (Fig. E]). About 3 cycles of wiggle are seen on each side 
of the source. The specific angular momentum for the orbital motion is ~ 3.7±1.3 AU km 
s~^, coni parable to that of the jet rotation yet to be confirmed, which is < 5 AU km s^^ 
( jPaper j ). 



Discussion 



4.1. SMMl: A very low-mass protobinary? 



It has been argued that almost all stars must form in bina ry or multiple systems because 
of the fragmentation of star-forming dense molecular cores (IGoodwin et al.ll2007| ). Could 
the source SMMl itself be a protobinary in the Class phase (i.e., the early phase of star 
formation) since the wiggle is reasonably reproduced by an orbiting source jet model? In 
order to investigate this possibility, we evaluate here the possible mass of the source, assuming 
that the source is a binary. Let Mj be the mass of the jet source, be the mass of the 
companion, and m be the mass ratio of the jet source to the companion, i.e., Mj = mMc, 
then the binary separation is a = (1 + m)Ro. From the Kepler's 3rd law of orbital motion, 
the total mass of the binary in terms of the model parameters is then 



9.5 X 10"^ (1 + mf 



K 



Az 



100 km s 



-1 . 



1°' 100 AU 



(7) 



With the best-fit values of n and Az, and the jet velocity Vj = 170 km s ^, we have 



6.9 (1 + Mj,p 



To determine the value of m, we resort to the previous estimation of the mass of the source. 
Previo usly, the mass of the s ource has been estimated to be ~ 60 Mj^p from an evolutionary 
model (IFroebrich et al.ll2003l ). This mass is also consistent with t hat deriv ed from the HCO^ 
rotating envelope-disk, assuming that the rotation is Keplerian (jPaper J l. If we adopt this 
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mass as the total mass of the binary, then m 1. This value of m is i n good agreement 



with that found in the binaries of very low- mass stars and brown dwarfs (IClose et al.l 12003 



Bouy et al .112 0061 ). which tend to have equal mass companions. In addition, the binary sepa- 



ration would be 2/?o ^ 4 .6 AU, also consistent with that found in those binaries (I Close et al. 



20031 : iBouy et al.l l2006l ). which have a range of separation between 0-30 AU. Therefore, 



SMMl itself could indeed be a very low-m ass protobinary becaus e of the fragmentation at 
the beginning of star formation (see, e.g., iMachida et al.l l2008bl ). It is unclear, however, 
why the jet source, with only ~ 30 Mj^p, can have such an energetic coUimated jet, but 
the companion as massive as the jet source does not have a jet or an outflow. It is also 
unclear if this protobinary will remain substellar in the final stage since it is surrounded 
by a compact envelo pe with an estimated mass of ~ 50 Mjup (altho ugh it is a lower limit. 



see 



Lee et al.l l2007bl ) at the eastern edge of the ammonia envelope (iPaper J )? Also, since 
the source SMMl could form a triple sy stem wit h the source SMM2 detected at ~ 84 AU 
to the southwest with a planetary mass (jPaper j ). could the source SMM2 be ejected from 
the SMMl binary sys tem in the early stage of star formation because of dynamical decay 
(iGoodwin et al.ll2007l )? The HH211, if indeed a triple system consisting of a close binary 
(<30 AU) with a t hird compo n ent a t ~ 100 AU, would be quite similar to those found in 
Upper Scorpius by iBouy et al.l (120061 ). 



4.2. The mass-loss rate and accretion rate 



Here we estimate the mass-loss rate of the jet and the accretion rate toward the central 
source, in order to compare their ratio with that predicted in the current jet launching 
models. The (two-sided) mass-loss rate of the jet can be given by 



(9) 



where N and h are the column density of the jet and the beam size perpendicular to the 
jet axis, respectively, as given in §2.3[ and is the mass of an H2 molecule. With the 



jet velocity Vj ~ 170 km s ^, the mass-loss rate is Mj ~ 1.8 x 10 ^ Mq yr ^, as found in 



Lee et al.l (l2007bl ). The mechanical luminosity of the jet is then 







(10) 



rough ly the same as the bolometric luminosity of the source, which is Lboi ~ 3.6L0 ( iFroebrich 
2005h. The jet is likely powered by a ccretion as in the current jet launching models (see, e.g., 
Shu et al.ll2000l : iPudritz et al.l 120071 ). and our result suggests that the jet can indeed carry 
away a large fraction of the accretion power in the early phase of star formation as pointed 
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out by ICabritI (120001 ) . Therefore, the accretion rate should be estimated using both the 
bolometric luminosity of the source and the mechanical luminosity of the jet. Assuming that 
the source SM Ml is a single protostar with a stell ar mass of ~ 60 M^„^ (IFroebrich et al. 
2OO3I : IPaper and a stellar radius of ~ 2Rq JstahlerlEgSsI : iMachida et aPboOSah . then 
the accretion rate is Ma ~ (i^boi + Lj)R^/GM^ ~ 8.5 x 10~^ Mq yr~^. Thus, the mass-loss 
rate is estimated to be ^ 20% of the accretion rate, con sistent with that predicted in the 



current MHD models fIShu et al.ll2000l : IPudritz et al.ll2007n . However, the source SMMl itself 
could be a binary (see §4.11) . If the jet source has a mass of ~ 30 Mjup and it is responsible for 
half of the bolometric luminosity, then the accretion rate toward the jet source is ~ 1.3 x 10^^ 
Mq jt~^. In this case, the mass-loss rate is ~ 14% of the accretion rate, also consistent with 



Paper 1 


(e.g., 


Soker 


2009) 



our estimation. 

Now let us check if the accretion rate is acceptable by comparing the accretion time with 
the outflow dynamical time. If the accretion rate was the same in the past, the accretion 
time would be M^/M^, i.e., ~ 7000 yrs if the source is single, or ~ 2300 yrs if the source is 
a binary. T he outflow dynamica l time can be estimated with the tip of the eastern outflow 
lobe in H2 (iHirano et al.l |2006| ) because the bright H2 emission there requir es the shock 
veloc i ty there to be < 55 k m s~^, otherwise H2 there will be all dissociated (jSmith et al. 
I991I : lO'Connell et al.l 120051 ). Therefore, with the tip at ~ 12500 AU from the source, the 



dynamical time for the outflow is > 1100 yrs. As a result, the accretion time is not shorter 
than the outflow dynamical time, and thus the accretion rate is acceptable. The accretion 
rate is high probably because of the extremely young age of the source SMMl. Note that, 
if the source SMMl is a binary, then in our estimation the accretion time is only two times 
the outflow dynamical time. If that is the case, the jet could be launched soon after the 
accretion started. 



4.3. Origin of SO and SiO emission: shock enhancement 

The abundances of SiO and SO in the gas phase are found to be highly enhanced in 
the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the 
source where the dynamical time scale is < 10 yrs. Here we discuss the possible origins of 
the abundance enhancement of these molecules. 

Although SiO abundance enhancement has been found in many Class jets, its origin 
is still not well determined. In HH 211, the abundance enhancement of SiO is closely related 
to the shocks, with the emission associated with a broad range of velocities. As pointed out 
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by Schilke et al. (1997) in their early work, the shocks are hkely required to generate the 
SiO gas by either (1) a release of SiO from dust grains, or (2) a release of Si atoms followed 
by conversion to SiO through gas-phase chemistry. The jet indeed seems to be d usty, wit h 
the 352 GHz continuum emissi on extending along the jet axis from the source (jPaper J ). 
Recently, iGusdorf et al.l ( l2008al ) have studied the formation time of the SiO from Si via gas- 
phase reactions in stationary C-type shocks with grain sputtering of Si-bearing material for 
the preshock number density of hydrogen nuclei nn = 10^ — 10^ cm~^. From their results, 
if the shock velocity is greater than 30 km s~^, we can approximate the SiO formation time 
to be roughly given by lO'^/nn yrs. If we extrapolate that SiO formation time for HH 211 
that has nn = 2n > 6 x 10^ cm~^ (see §2.31 for n), then the formation time is < 2 yrs 
for HH 211, short enough to produce the SiO emission in the SiO knots. Note that non- 
stationary C-ty pe shocks ([G usdor f et al.ll2008bl ) and J- type shocks with more sophisticated 
grain processes (iGuillet et al.ii2009l) could also be responsible for the SiO emission and more 
studies are needed to explore those possibilities. There is a lack of SiO emission [not only 
i n J=8-7 transition, but a lso in lower transitions in J=5-4 (Gueth et al. in prep) and J=l-0 
( IChandler fc Richerll200l[ )] at the base of the jet in the innermost pair of CO knots, BKO 
and RKO. This lack of SiO emission is unlikely due to a lack of material because the jet is 
expected to have the same density along th e jet axis except near the launching point where 
the density is much higher ( IShu et al.ll2000l ). This lack of SiO emission suggests that the SiO 
is not abundant there in the gas phase. If the SiO were there abundant in the gas phase, 
the shocks tha t excite the CO emis sion (which was also detected in higher transition in the 
far infrared by lGiannini et al.ll200l[ ) would also excite the SiO emission. 



SO abundance enhancement has also been seen toward other Class sources, e.g., L1448, 
HH 212, and L1157, because of shocks. In L1157, the SO emission is seen in the prominent 
bow shocks at > 20000 AU away from the source, but no jetlike SO emission is seen around 
the source ( Bachiller et al.l 2001). In L1448, the SO emission is seen at shock-precursor and 
post-shock velocities ( iJimenez-Serra et al.ll2005l ). In HH 212, the S O emission at the base is 
jetlike as in HH 211 and it extends to ~ 800 AU from the source (iLee et al.ll2007a] ). Note 
that in that source SO emi ssion in lower t ransition has also been detected upto 8000 AU 
to the south along the jet (ILee et al.ll2006l ). Thus, the abundance of SO in the jet seems 
to decrease from the younger sources, e.g., HH 211 and HH 212, to the older sources, e.g, 
LI 157 (which shows no clear jetlike SO emission), like that of SiO (jShang et al.l 120071 ). In 
LI 157, it was proposed that sulfur is rel eased from dust grai ns in the form of H2S in the 
shocks and that is then oxidized to SO ( iBachiller et al.ll200ll ). In L1 448, H2S emission is 
presen t in the shock-precursor, but is missing in the post-shock gas. iJimenez-Serra et al. 
( I2OO5I ) consid ered oxidation o f H2S to SO in shock, but judged the timescale [100 - 1000 yrs 
suggested by IWakelam et al.l (|2004[ )] too long. They thus proposed that the SO molecules 
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might be abundant on dust grains and directly released from grains in the shocks. HH 211 
and HH 212 could have the same origin of abundance enhancement in SO. In HH 211, the 
shocks can release sulfur from grains in the form of e.g., H2S, and then oxidize it to SO, as 
proposed in L1157. If the je t is rich in atornic ox ygen (01) and carbon ion (CII) (as can be 
seen in ISO observations in iGiannini et al.ll200ll ). then oxidation of H2S to SO could take 
only a few years. For example, the following sequence would have timescales less than 5 
yrs in a shock environment with nn ~ 10^ cm~'^, Tgas ~ 1000 K, n(HI)/n(H2) ~ 0.01, and 
n(OI)/nH ~ 10-^: H2S + H ^ HS + H2; HS + O ^ SO + H. In this exa mple, the rate 



coefficients are 2.5x10 cm s 



^ and 1.3x10"^° cm'^ s"^ 



(RATE06 by IWoodall et al. 



20071 ). and the conversion timescales for H2S — ?■ HS and HS — t- SO are 0.13 yrs and 3 yrs, 
respectively. (The timescale for HS + H S + H2 is 0.13 yrs, but the timescale for reverse 
S + H2 — > HS + H is 2.6 yrs at 1000 K. Thus, about 5 percent of sulfur atoms will be in HS 
in 3 yrs and be available for HS + O.) Again, there is a lack of SO emission at the base of 
the jet in the innermost pair of CO knots. This suggests that the SO is not abundant there, 
neither in the gas phase nor on the dust grains. If the SO were there on the dust grains, it 
would evaporate into the gas phase and then be excited to produce the SO emission. 

We now speculate what could cause the lack of SiO and SO abundances in the gas phase 
at the base of the jet in the innermost pair of CO knots. Since the mass-loss rate in the jet is 
high (~ 1.8 X 10~® Mq yr~^), SiO and CO as well could have form ed via gas-phase react ions 
in an initially atomic jet close to the launching point (< 0.1 AU) ( IGlassgold et al.lll99lh . If 
that is the case, the lack of gas-phase SiO in the innermost pair of CO knots suggests that 
the SiO could be afterward dissociated e.g., by the possible presence of a far-UV radiation 
field near the source. Howe ver, it is also possible that the SiO is rapidly depleted onto the 
grains or converted to Si02 (iGusdorf et al.ll2008al ) if the temperatur e of the jet indeed drop s 
rapidly with the distance from the launching point as predicted by IGlassgold et al.l ( 1l99ll ). 
In addition, the lack of gas-phase SO in the innermost pair of CO knots also suggests that 
the SO, if formed at the base, is afterward destroyed. On the other hand, the CO could 
stay in the gas phase, giving rise to the CO emission even in the innermost pair of CO 
knots. Both the SiO and SO in the gas phase can form again in the shocks, as discussed 
above, producing the SiO and SO in the downstream. The SO and SiO abundances in the 
jet are seen decreasing with the outflow dynamical age, probably because the mass-loss rate 
decreases with time and so does the formation rate of the SiO and SO molecules. 
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5. Conclusions 

SiO, CO, and SO can all be a good tracer of the HH 211 jet, tracing the internal shocks 
in the jet. Although the SiO, CO, and SO emissions show roughly the same morphology of 
the jet, there are detailed differences. In particular, the CO emission traces the jet closer 
to the source than the SiO and SO emissions. In addition, in the better resolved internal 
shocks, the CO emission is seen slightly ahead of the SiO emission, likely because the SiO 
emission traces stronger shocks closer to the shock front and the CO emission traces the 
weaker shocks in the downstream. In the internal shock where the SO emission is bright, 
the SO emission appears slightly ahead of the SiO emission, tracing weaker shocks than the 
SiO emission. 

The jet is seen with more than one cycle of wiggle on both sides of the source SMMl in 
SiO, CO, and SO. The wiggle is reflection-symmetric about the source and can be reasonably 
fitted by an orbiting source jet model. The best-fit parameters suggest that the source itself 
could be a very low-mass protobinary with a total mass of ~ 60 Mjup. The protobinary may 
have two equal-mass sources and a binary separation of ~ 4.6 AU, in good agreement with 
those found in the binaries of very low-mass stars and brown dwarfs. 

The abundances of SiO and SO in the gas phase are found to be highly enhanced in 
the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the 
source where the dynamical time scale is < 10 yrs. The abundance enhancements of these 
molecules are closely related to the internal shocks in the jet. The detected SiO is either the 
consequence of the release of Si-bearing material from dust grains or of its formation via gas 
chemistry in the shocks. The SO, on the other hand, seems to form via gas chemistry in the 
shocks. 

We thank the SMA staff for their efforts in running and maintaining the array, and 
the anonymous referee for the precious comments. C.-F. Lee thanks Prank Shu for fruitful 
conversations and Ngoc Phan-Bao for useful discussion. A. P. is grateful to Joaquin Santiago- 
Garcia, Robert Estalella and Nuria Huelamo for useful discussions. 
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Fig. 1. — SO and CO maps of the jet on top of the SiO map (gray cont ours) of the jet and 
the 352 GHz continuum map (green contours) of the envelope-disk from iPaper t . The maps 
have been rotated by 63.4° counterclockwise. The blueshifted and redshifted SiO emission 
are integrated from -21.2 to 9.2 km s~^ and from 9.2 to 47.5 km s~^, respectively. The cross 



marks the position of the source SMMl. The knots have the same notations as in iPaper 1 



(a) shows the blueshifted and redshifted SO emission, integrated from -16 to 9.2 km s~^ and 
from 14.5 to 43.5 km s~^, respectively. No SO emission is detected from 9.2 to 14.5 km/s. 
The contours start at 0.52 Jy beam~^ km s^^ with a step of 0.52 Jy beam~^ km s~^. The 
SiO contours start at 1 Jy beam~^ km s~^ with a step of 1.5 Jy beam~^ km s^^. (b) and 
(c) show the blueshifted and redshifted CO emission, integrated from -16 to km s^^ and 
from 20 to 43 km s~^, respectively. In (b), the CO contours start at 0.56 Jy beam^^ km 
s~^ with a step of 0.84 Jy beam~^ km s~^. The SiO contours are the same as in (a). In 
(c), the CO contours start at 0.5 Jy beam~^ km s~^ with a step of 0.5 Jy beam~^ km s~^. 
The SiO contours start at 0.84 Jy beam~^ km s~^ with a step of 0.84 Jy beam~^ km s^^. 
The synthesized beams for the lines are 0''46x0''36 (natural weighting) in (a) and (b), and 
0!'35x0''25 (super-uniform weighting) in (c). 
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Fig. 2. — PV diagram cut along the jet axis in CO (red image and contours) for the two 
inner pairs of knots, e.g., B KO, BKl , RKO and RKl, in comparison to that in SiO (green 
image and contours). As in iPaper I . the yellow lines mark the positions of the sub-knots 
in knots BKl and RKl, and the question marks indicate the possible positions of the two 
sub-knots in knot BKl. The blue lines mark the positions of knots BKO and RKO. The 
systemic velocity is 9.2 km s~^. The white horizontal dashed lines mark the mean velocities 
of the jet, which are —6 km s~^ on the blueshifted side and 28 km s~^ on the redshifted 
side, as given in iPaper 11 . The SiO contours start at 0.15 Jy beam~^ with a step of 0.15 Jy 
beam~^. The CO contours start at 0.1 Jy beam~^ with a step of 0.1 Jy beam~^. 
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Fig. 3. — PV diagrams cut along the jet axis in CO (red image and contours), SiO (green 
image and contours), and SO (blue image and contours) for knots BK2 and BK3. The white 
dashed line indicates the mean velocity for the eastern component of th e jet. T he orange 



dashed lines separate the heads and tails for knots BK2 and BK3, as in iPaper J . The SiO 
contours start at 0.15 Jy beam~^ with a step of 0.15 Jy beam"^. The CO contours start at 
0.1 Jy beam~^ with a step of 0.1 Jy beam~^. The SO contours start at 0.08 Jy beam~^ with 
a step of 0.08 Jy beam~^. 
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Fig. 4. — Schematic diagram of the orbiting source jet model. Star marks the source position 
in the orbit. 
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Fig. 5. — The orbiting source jet model for the wiggle plotted on top of Figured] The red 
and blue curves are derived from our best-fit model (see the text in 



